THE NORMAL BEHAVIOR OF 
THE SMARANDACHE FUNCTION 


KEVIN FORD 


Let S(n) be the smallest integer k so that n|k!. This is known as the Smarandache 
function and has been studied by many authors. If P(n) denotes the largest prime 
factor of n, it is clear that S(n) > P(n). In fact, S(n) = P(n) for most n, as noted 
by Erdés [E]. This means that the number, N(z), of n < x for which S(n) 4 P(n) 
is o(z). In this note we prove an asymptotic formula for N(z). 

First, denote by p(u) the Dickman function, defined by 


* p(v—1) 


du ys 
F v (u>1) 


p(u)=1 (0<u< ll), pu) =1- f 


For u > 1 let € = €(u) be defined by 
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It can be easily shown that 


logs u 
€(u) = logu+loggu+O ( ee ) 


where log, z denotes the kth iterate of the logarithm function. Finally, let ug = 
uo(z) be defined by the equation 


log x = up€(uo). 
The function uo(z) may also be defined directly by 
log x = uo Cue — 1) : 
It is straightforward to show that 
a) wo= (BEE) (1 He + 5B +0 ((BBE) J). 


We can now state our main result. 
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Theorem 1. We have 


N(a) ~ Vi(1 + log 2) 


93/4 (logz log, x)9/4z} 1/40 (ug). 


There is no way to write the asymptotic formula in terms of “simple” functions, 
but we can get a rough approximation. 


Corollary 2. We have 


N(x) = zexp {-(v2+ o(1))./log x logs z} . 


The asymptotic formula can be made a bit simpler, without reference to the 
function p as follows. 


Corollary 3. We have 


log z 
7(1 + log2 ne eM —1 
N(e) ~ SUPE) tog 2) (logs 2) 2/40 exp [ 6 otal, 


where y = 0.5772... is the Euler-Mascheroni constant. 


This will follow from Theorem 1 using the formula in Lemma 2 which relates 
p(u) and €(u). 

The distribution of S(n) is very closely related to the distribution of the func- 
tion P(n). We begin with some standard estimates of the function U(z,y), which 
denotes the number of integers n < x with P(n) < y. 


Lemma 1 {HT, Theorem 1.1]. For every « > 0, 


log(u + 1) ss log x 
log y ~ logy’ 


U(z,y) = zp(u) ¢ +0 ( 


uniformly in1 <u < exp{(logy)9/> ¢}. 
Lemma 2 [HT, Theorem 2.1]. For u > 1, 


(-0(2)) Harb feoa} 
= {-u (1084+ og u 21H Ga) p 


Lemma 3 [HT, Corollary 2.4]. [fu> 2, |v] < u/2, then 


p(u) 


p(u—v) = p(u) exp{v€(u) + O((1 + v?)/u)}. 
Further, ifu>1land0<vu <u then 


p(u—v) « p(ues™. 
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We will show that most of the numbers counted in N(z) have 


P(n) ~ exp { \/Flog tog, =} : 


Vi = ep {}Vbogzlons}, Ya=¥i = exp {2Viogz Toes}. 


Let N, be the number of n counted by N(z) with P(n) < Yi, let No be the number 
of n with P(n) > Yo, and let N3 = N(x) — Ni — No. By Lemmas 1 and 2, 


Let 


N, < U(z, Yi) = rexp{—(1.5 + o(1)) log z log, zr}. 


For the remaining n < xz counted by N(z), let p = P(n). Then either p?|n or for 
some prime q < p and b > 2 we have q° || n, q’ { p!. Since p! is divisible by q'?/¢! 
and b < 2logz, it follows that g > p/(3logz) > p'/?. In all cases n is divisible by 
the square of a prime > Y2/(3logz) and therefore 


6z log z 
Nex S- as Y « cexp {-1.9\/logz log 2}. 


¥2 
P2 Flog 


Since q > p'/? it follows that q'?/4 || p!. If n is counted by N3, there is a number 


b > 2 and prime g € [p/b, p} so that g°|n. For each b > 2, let N3 (x) be the number 
of n counted in N3 such that q° || n for some prime g > p/b. We have 


31 : 
S_ Ns. <x (*) < nexp {-(5/3 + 0(1)) Vlog zlog, =} . 
1 
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Next, using Lemma 1 and the fact that p is decreasing, for 3 < b < 5 we have 


sen 5, ((59)+ E964) 


Yi <p<Y2 p/bxq<p 
1 log z 1 log z — logp — blogg 
«zc \- oP (Zs - ) +» ag? (epee nee ieee 
Yi<p<Y2 2 p/2<q<p = 
log x 
«zt 5° p °o (2 -(b+1)), 
Yi<p<Yo SP 


By partial summation, the Prime Number Theorem, Lemma 2 and some algebra, 


N3,5 < exp {-(15 + 0(1))/log z log, x} : 
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The bulk of the contribution to N(z) will come from N39. Using Lemma 1 we 
obtain 


(2) 
won, oS) ol) 


Y, <p<Yo 5 <a<p 
p (282 ) p (282 lege _ 2) 
" logs z log p log g 
=(1+0(, a) )2 > p2 - 2 pq? 
1<p<Y2 p/2<q<p 


By Lemma 3, we can write 


gxr— yt 
p (WES PEP 2) = p (PEE 3) (1+0[ 22) | 
log q log g log z 
The contribution in (2) from p near Y; or Y2 is negligible by previous analysis, and 
for fixed q € [¥1, Yo/2] the Prime Number Theorem implies 


Y= PE + oliogs) *) = $2 +0 (3), 


2 
Rope ase log g log p log* Y; 


Reversing the roles of p,q in the second sum in (2), we obtain 


1 log x log2 /logz 
ry = logs z as _9 — _ ; 
vn € uy ( logs )) + 2 p (> (Be ) * jogp” (Pee 3)) 


Yi<p<Yo2 


By partial summation, the Prime Number Theorem with error term, and the change 
of variable u = logz/logp, 


B. Ns= (1 +0 (#2) a (a) ii elu -3))a 1/u gy 


where 
logz 


U=r uz = bu}. 


log, x’ 


The integrand attains its maximum value near u = uo and we next show that the 
most of the contribution of the integral comes from u close to up. Let 


log; z ae 
log, z 


iit w,= Kuo, we=wi ( 


where K is a large absolute constant. Let J; be the contribution to the integral in 
(3) with |u — uo| > wo, let Iz be the contribution from w; < |u — ug| < wo, let J3 
be the contribution from w2 < |u — up| < wi, and let J, be the contribution from 
|u — uo| < we. First, by Lemma 2, the integrand in (3) is 


wo{-(L-feau) vgs}. (8) 
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The function 1/c + c/2 has a minimum of V2 at c = V2, so it follows that 


re exp {- (v2+10 *) Vlogs log, =}. 


Let u = up — v. For w; < |v| < wo, Lemma 2 and the definition (1) of up imply 
that the integrand in (3) is 


log z vv xs v? 
< — Sn pas as 
< p(uo) exp {r6(w) € ahora tg 5] +O (= + log uo 


v? vu? 
< plu) 1/*° exp {- ages +O (= + log )} 
0 


2 
< p(uo)z /*¥° exp {70.955 log x} 
0 
for K large enough. It follows that 
In < uop(uo)z /%° exp{—20 log, z} « (logz) ?°p(uo)z Tf Mo: 
For the remaining wu, we first apply Lemma 3 with v = 2 and v = 3 to obtain 


Untwy 2£(u) log 2 
I3+Ig4={1+0O / 282z ) p(ujx V4 : 4+ 28% 364) ) du 
Ogt a u logz 


0 Wi 


We will show that J3 + I, >> p(uo)z 1/9 (log z)3/?, which implies 


uoptwi 2€(u) 
(4) N(z)= (1+0/ ) i: p(u)z mis a sbe eet) du. 


abs u * loge 


This provides an asymptotic formula for N(x), but we can simplify the expression 
somewhat at the expense of weakening the error term. First, we use the formula 


= logs u 
E(u) = logu + loggu+O (72) : 


and then use u = uo + O(u Ae and (1) to obtain 


uoptw, 


Igt+I,= (a +O (les) 2 (1 + log 2)z(log z)? (log gz)? / p(u)z 1/u Ay. 


QO Wi 


By Lemma 3, when wz < |u| < wi, where u = up — v, we have 


yee 1 ] 2 a) 
plug -v)z 0 * <€ p(ug)z *o exp < vé(ug) — ee (= Cae =) 


As 
<K p(uo)r “0 exp 


< p(uo)z *o exp {-5 log z} 


« p(uo)z *o(logyz) 3 


provided K is large enough. This gives 
j pluje /* du < p(uo)e ™*°(logz)!/4(logy x) 35. 
w2K|u uo[<wi 


For the remaining v, Lemma 3 gives 


plug —v)z Vo ») — (a +O (3) p(uo)x 1/¥° exp {-S tog} 


Therefore, 
uptw2 we 21 . 
p(uo) 1pue / p(ujz V%du= (a +O (22) / exp {-9 ae} dv. 
ug we we 0 


The extension of the limits of integration to (—oo,0o) introduces another factor 
1+ O((logz) 1), so we obtain 


seve m(1 + log 2) tp 
fat te= (140 (BEE)) are Coes ton, 2)"p(uole 


and Theorem 1 follows. Corollary 2 follows immediately from Theorem 1 and (1). 
To obtain Corollary 3, we first observe that €’(u) ~ u | and next use Lemma 2 to 


write 
a i: ” €(t) dt 
p(uo) Jiri €xp ; E( . 
By the definitions of € and ug we then obtain 
uo (uo) aan 
/ ewar= f wise ro 
1 0 
co 
= ef(uo) _ vb 
) 


log =z 
log z =i a is! 


Uo 


Corollary 3 now follows from (1). 
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